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5^ ■ Abstract 



The ground state energy and the extend of the wavefunction of a negatively charged donor {D ) 
located near a semiconductor-metal or a semiconductor-dielectric interface is obtained. We apply 



'^ ' the effective mass approximation and use a variational two-electron wavefunction that takes into 

^ ■ account the influence of all image charges that arise due to the presence of the interface, as well as 

the correlation between the two electrons bound to the donor. For a semiconductor- metal interface, 
K* ' the -D^ binding energy is enhanced for donor positions d > l.bas {as is the effective Bohr radius) 

in ■ 



due to the additional attraction of the electrons with their images. When the donor approaches the 
interface (i.e. d < l.Sa^) the D~ binding energy drops and eventually it becomes unbound. For a 



in 

^D I semiconductor-dielectric (or a semiconductor- vacuum) interface the D binding energy is reduced 



for any donor position as compared to the bulk case and the system becomes rapidly unbound 
when the donor approaches the interface. 
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I. INTRODUCTION 

The negatively charged donor center, also called D~ system (a hydrogenic donor with a 
second electron bound to it) has received a lot of interest in the past both from researchers 

n n 

in astrophysics (where it is known as the H ion [1]) and semiconductor physics [2|. Further- 
more, the D~ center is one of the simplest "many-body" electronic system that can also be 
used as a model system to test how well certain theories are able to include electron-electron 
(e-e) correlations. It has been shown that in the absence of electric and magnetic fields the 
D^ has only one bound state |2|]. Experimentally, D~ states have been observed in bulk 
semiconductors |3|, in quantum wells j^, and in superlattices |5[. Tunneling through the 
D^ state was reported in Ref. la for a D~ confined in a double-barrier resonant tunneling 
device. The D~ resonance appears in addition to the known resonance due to tunneling 
through the ground state of the neutral donor D^. It was found that in high magnetic fields 
the amplitude of the D~ resonant peak becomes significantly larger as compared to the 
amplitude of the Z)° peak [6|. 

Theoretically, the D~ center has been investigated using the variational method within 
the effective mass approach: i) for the bulk case 7[], ii) the two-dimensional case [sl, |9[, 
and iii) for the superlattice case [10|. In a double-quantum- well device a neutral donor 
D^ at the center of the quantum well in the presence of magnetic field can bind a second 



electron in the other well (a spatially separated D~ center) as shown in Ref. [ll|. In the 
case of the remote D~ center, where donor and electrons are located in different quantum 
wells, strong electron-electron correlations can give rise to magnetic-field-induced angular- 



momentum transitions 
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12| . It has been predicted that electron-phonon interaction in a 



weakly polar semiconductor leads to a substantial increase of the D~ binding energy 13] . 
The negatively charged donor center was used as an approximate model system for a trion 
(charged exciton) to explain the experimental behavior of the two-dimensional electron gas 
in a quantum well in high magnetic field at high- laser power [141 ]. 

In the last years, there has been renewed interest in the study of dopants in semicon- 
ductors due to the possibility to dope the material in a controlled way and to tailor the 
electronic properties in order to create new functional devices. Because of the increased 



miniaturization, the dopant atoms appear closer and closer to interfaces 15|. Recently, us- 



ing STM, the binding energy of individual dopants close to a semiconductor interface was 



measured and found to be substantially increased [15J]. In a recent transport experiment on 
a nanowire surrounded by a metallic gate it was suggested 16| that signatures of the D~ 
state were observed. Due to the closeness of the metallic gate to the donor it was argued that 
the metallic gate screens the repulsive e-e interaction which should lead to a larger binding 
energy. Here we will show that this expectation is only correct if the donor is not too close to 
the interface. For very close proximity of the donor to the interface the D~ binding energy 
drops and can even become negative, i.e. unbound D" system. A second motivation for the 
present study is that the D~ system can serve as an entangled pair of electrons which cur- 



rently is of great interest for quantum information applications |17|. The quantum control 
of a D^ near a semiconductor-dielectric interface and its possible application for quantum 
computing has been discussed in Ref [l8|. 

In this paper we study the spin-singlet state of a D^ system near a semiconductor-metal 
(and semiconductor-dielectric) interface within the effective mass approach. As compared 
to the 3D situation the present problem differs in the following two aspects: 1) the many- 
particle wavefunction is zero at the interface, and 2) due to the dielectric mismatch at the 
interface image charges are induced that results in a complicated multi-center Coulomb 
problem. For the ground state trial function we used a Chandrasekhar type space sym- 
metric wave function which we modified in order to satisfy the boundary condition on the 
semiconductor-insulator interface, and to take into account the contribution of the Coulomb 
interaction with the image charges in the system. To describe the interaction of the elec- 
trons with the images in the insulator or in the metallic gate, we add additional terms to the 
Chandrasekhar trial function [l|]. For the case of a single electron bound to a donor near a 
metallic or dielectric interface, i.e. the D^ problem, similar terms were introduced recently 



in Ref. 
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which we modified slightly in order to obtain even better results for d> as- In the 
present paper, we have calculated the ground state energy of the D~ system as a function 
of the donor position with respect to a semiconductor-metal, a semiconductor-dielectric and 
a semiconductor-vacuum interface. 

The present paper is organized as follows. In Sec. Hi] we present the two electron Hamil- 
tonian and propose trial wave functions for the D" center. In Sec. IIIII the ground state 
energy of the D~ center near different interfaces is studied as a function of the position of 
the donor. The extend of the wavefunction and its average position in the direction parallel 
and perpendicular to the interface were also calculated. Our conclusions and a summary of 



our results are presented in Sec. [IVl 

II. THE FORMALISM AND THE VARIATIONAL WAVE FUNCTION 

The Hamiltonian of the D~ system, consisting of a donor at a position r^ = (0, 0, d) near a 
semiconductor-metal (semiconductor-dielectric) interface and two electrons is, in cylindrical 
coordinates, given by the expression 

H = H^ + H, + U{r^,r2), (1) 

where 



2 



d^ Id 1 d^ d^ 

dpi Pi dp, pfdef dz^ 



I J 



+ f/°(?.) (2) 



is the Hamiltonian of a neutral D center near an interface, with 

f/0(f,) = ^ - ^ ^ - ^ ^ (3) 

the Coulomb interaction terms and where dimensionless units of the effective Bohr radius 
as = f?e\jm,^e^ and twice the effective Rydberg energy 2i?* = h^ lm±a?g were used and 
o" = mx/mw is the ratio between the transverse and longitudinal effective mass (in our 
numerical calculations we assume a = 1). 

The first term in f/°(rj) describes the interaction between an electron and its image, the 
second arises due to the attractive interaction between an electron and the donor, and the 
third term is due to the interaction between an electron and the donor image (as well as 
the donor and the electron image). See Fig. [T] for a schematic representation of the different 
Coulomb terms in case of a semiconductor-metal interface for which Q = —1 (i.e. Eq — ^ oo). 
In Eq. ([3]) the image charge is given hj Q = {es — £o)/{^s + £o) with es {£o) the permittivity 
of the semiconductor (dielectric). For the case of a semiconductor- metal interface we assume 
a very thin oxide layer between the semiconductor and the metal and its only effect is to 
prevent the electron to penetrate into the metal, i.e. it provides a very high potential barrier. 
The two-electron Coulomb potential has the following form 

f/(ri. Fa) = f/ee(ri. Fa) + Uei{ri, Fa), (4) 



with 

f/ee(ri,?2) = — , (5a) 

a/p? + pi - 2pip2 cos(6'i - 6*2) + {zi - Z2Y 

f/ej(?l,?2) = — , (5b) 

Vp? + pi - 2pip2 cos(^i - Q2) + {zx + ^2)^ 

where f/ee describes the electron-electron interaction and Vei is the interaction between an 
electron and the image of the other electron. 

The potential energy for the D° electron along the z-direction for p = is shown by 
the dashed curve in Fig. [2]^a) for the semiconductor-metal interface and in Fig. EJ^b) for 
the semiconductor-dielectric interface. The electron is subject to a — e/|r — rd| Coulomb 
potential near the donor and the potential due to the image charges. When an electron is 
bound to the positive donor the system becomes neutral and this electron will screen the 
donor Coulomb potential for a second electron. The mean-field potential seen by the second 
electron is given by 

t^°(?2) + (^^o(ri)|f/(ri,r2)|^z3o(?i)) (6) 

where ?/'£)o(r) is the normalized ground-state wavefunction of the D° electron which we 
obtained previously in Ref. il9|. This potential is plotted in Figs. [2](a,b) by the solid curve 
and it is immediately clear that the second electron will feel a strongly screened donor 
potential. The 1/r Coulomb potential is replaced 8| by the screened potential (l + r)e~^''"/r. 
As a result the second electron can at most only be very weakly bound. Solving this mean- 
field problem for the second electron will strongly underestimate the D~ binding energy 
because it neglects: i) exchange, and ii) electron-electron correlation. The latter will lead 
to a polarization of the neutral D^ system. 

In order to account for exchange and correlation we introduce a Chandrasekhar type space 
symmetric wave function (because we consider the spin single state, i.e. V'(n, "^2) = "^(^2, "H,)) 
which also takes into account the interaction with the image charges in the system, and 
satisfies the boundary condition at the interface: 

V^(ri, r2) = N'p{zx) p{z2) [fi'(ri, r2) + g{r2, ri)]v?(ri, r2), (7) 



where 

p{zi) = Zi/{l + azi), (8a) 

5f(ri, r2) = exp(-Airi - \2r2) exp(-/3ip(zi) - f32p{z2)), (8b) 

V9(ri,r2) = l + (5ri2, (8c) 

with Tj = a/p^ + {zi — dy the distance between the electron and the donor, ri2 = 
■\/Pi + P2 ^ 2pip2 cos(^i — 6^2) + {zi — z^)'^ is the distance between the two electrons and 
A^ is a normalization constant. In our previous work [l9(] on the neutral donor problem 



we found that introducing an asymmetry in r,, i.e. Vi = \/ p\ + 7(-2j — (i)^, with 7 an ex- 
tra variational parameter, improved the binding energy only with 2.7% for dja^ = 1 and 
with less than 1% for d/as = 5. Therefore, in order to limit the number of variational 
parameters we took 7 = 1. The factor p{zi) in front of the wave function ([7]) guarantees 
that the boundary condition ip{zi = 0) = is satisfied for each electron and that for large 
d-valnes the 3D result can be recovered (because then p{zi) approach the constant value 
1/a). The wave function contains 6 variational parameters: Ai, A2, /3i, (32, S and a. The 
symmetric combination of hydrogen-like factors exp(— Ajr^) in Eq. (18bp . takes into account 
the interactions of the two electrons with the donor, as was proposed in Ref. qJ, and V5(ri, r2) 
describes the electron-electron correlation. We introduced the combination of exponential 
factors f{zi) = exp{—l3ip{zi)) that describes the overall interaction of each electron with its 
image, as well as with the images of the donor and the other electron. Similar functions for 
the electron-electron repu 



system in quantum wells 



sive interaction (but with a = 0) were used previously for the D~ 
lOJ I and for the two-electron parabolic quantum dot 20 1. 
For a neutral electron-donor D^ problem near an interface 19[ we took a = and thus 
p{z) = z which resulted in very good agreement between the results of the variational 
method and a numerical 'exact' finite-element solution. Although the energy was found 
to be accurate it is clear that the bulk wavefunction is not recovered ioi d ^ 00, which 
should be spherical symmetric around the donor. This is the reason why we introduced 
p{zi) = Zi/{1 + azi) which for 2; — )■ 00 approaches the constant value 1/a and makes the 
wavefunction spherical symmetric. In Table H] we compare the obtained variational energy 
for the D~ using p{zi) = Zi and p{zi) = Zi/{1 + azi); we see that the latter gives lower 
energy for all (i-values. Thus this extra variational parameter has only a small influence on 
the energy (we checked that the same conclusion holds for the D^ problem). But we found 



d/aB 


P{Zi) = Zi 


p{z 


E/2R* 
i) = Zi/{1 + azi) 


asymmetric e-e 


2 


-0.40883 




-0.41464 


-0.41464 


3 


-0.49106 




-0.49783 


-0.49783 


4 


-0.51680 




-0.52180 


-0.52186 


5 


-0.52518 




-0.52885 


-0.52891 


6 


-0.52768 




-0.53061 


-0.53067 


7 


-0.52811 




-0.53059 


-0.53063 


8 


-0.52787 




-0.52998 


-0.53001 


9 


-0.52747 




-0.52922 


-0.52924 


10 


-0.52709 




-0.52850 


-0.52852 



TABLE I: The ground state energy of a D~ center near a semiconductor-metal interface at different 
donor-interface distance d, calculated by using two different p{zi) terms in the wave function, 
Eq. d?]). In the last column we present the results for a 7^ and including the asymmetric e-e 
correlation term, Eq. Q, in the variational wave function. 



that quantities as, e.g., the average electron position behave much better at large d if we 
include a as a variational parameter. For rf — t- 00 we have p — )■ 1/a and our variational 
wavefunction reduces to the one proposed by Chandrasekar [l| which has been shown to 
result in accurate values for the bulk D~ binding energy. 

Due to the presence of the interface the e-e correlation term is not necessarily circular 
symmetric. In order to check the influence of this asymmetry on the energy we introduced 
an additional variational parameter 6z in the e-e correlation part of Eq. ([7]), 



ipniri, rs) = 1 + 6\ pj + pj - 2pip2 cos(6'i - 62) + S;,{zi - Z2Y. 



(9) 



The ground state energy of the D~ near a semiconductor-metal interface calculated by 
using (y9£) in the wavefunction Eq. (jTj) is presented in the last column of Table [B Notice that 
allowing the e-e correlation to be anisotropic does not have a significant influence on the 
energy, the differences between the energy calculated with or without 6^, are very small (less 
than 0.02%). Therefore, in the following we did not include 5^ in our calculations. 



The D binding energy is defined as follows: 

Eh{d) = Eno-En-, (10) 



where E^o is the ground state energy of the D^ system 19|. Eq. ( ITOj) gives the energy 
that is needed to remove one of the electrons from the donor to infinity. This definition 
corresponds to the one used for a 3D D^ system. For a metallic interface, another definition 
of the binding energy can be used: 

Ef,{d) = Ef + EDO-ED~, (11) 

where Ef/2R* = —0.0312 is the binding energy of an electron with its image near a metallic 



interface 
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2l|. This is the energy needed to remove one electron from the D~ and bring 
it to the interface. 

The expectation values for the position of electron i in the z-direction and in the p plane 
are calculated as (pj) = f Pi X12 dr and {zi — d) = f {zi — d) Xu dr, respectively, using only 
the first part in the total trial function, i.e. X12 = N p{zi)p{z2) 5'(ri, r2)v9(ri, r2). This allows 
us to discriminate between the two electrons, one is closely bound to the positive impurity 
while the other one is very weakly bound. Notice that thanks to our variational wavefunction 
we are able to separate 'artificially' the two electrons. But the variational wavefunction itself 
considers the two electrons to be indistinguishable because ip{^i,^2) = V'(?2,?i)- 

III. THE GROUND STATE ENERGY OF THE D CENTER 

A. Semiconductor-metal interface 

The ground state energy, in units of 2R*, is shown in Fig. [3] as a function of the distance 
of the donor from the interface d/aB- The D~ energy exhibits a shallow minimum around 
d ~ 6.5aB and approaches the bulk result Ebuikf^R* = —0.5259 from below for d — )■ cxd. For 
d < 6.5aB the D~ energy is a decreasing function of d which is mostly a consequence of the 
extra constraint that the wavefunction has to be zero at the interface. This can be inferred 
from the result without image charges (dashed green curve in Fig. |3]). The outer electron 
of the D~ center becomes bound for d > l.la^. We considered the second electron bound 
if the variational parameter A2 is larger than zero and the average position of the electron 
is not too far from the donor, i.e. {z2 — c?) ^ lOa^. For a semiconductor-metal interface, 

8 



we have two positive image charges from the electrons and only one negative image charge 
from the donor, so a strengthened electron bound state is expected. This is confirmed by 
the fact that if we ignore the image charges the second electron becomes only bound when 
d > 3.2aB (see the dashed green curve in Fig. [3]). 

The binding energy is plotted in Fig. H] using the two different definitions of the binding 
energy (Eqs. (fTO!) and (TTTl) ). A clear maximum is found for d ^ S.Sas beyond which the 
binding energy slowly decreases to its (i — )■ oo value. This algebraic decrease is due to the 
interaction with the image charges and is largely a consequence of the d-dependence of the 
D° energy. For the definition of the binding energy with Ef = (blue dashed curve in 
Fig. Hj) the bulk D~ binding energy Ehhuikf^R* = 0.0259 is reached from above. Notice that 
for d = 3.8aB, Ei,/2R* = 0.0823 which is a factor 3.18 larger than the bulk D~ binding 
energy. The large d-range is shown in the inset and we found that it can be fitted by the 
curve Eb/2R* = 0.0259 + 0.2550^/^. 

It is remarkable that in the D~ case the effect of the image charges on its energy almost 
cancels out for intermediate d, while this is not so for the neutral D^ system. In our previous 



work [19|] we found that in the case of a semiconductor-metal interface the contribution of 
the image terms to the energy of a neutral D^ center is given by AE = 1/Ad for large 
d. This is the sum of the contributions arising due to the interaction of the electron with 
its image (about —l/4d) and the interactions of the electron with the donor image as well 
as the electron image with the donor (each one is about 1/Ad). In Fig. [1] we show by 
dashed lines the repulsive interactions between the electrons and the donor image as well as 
between the donor and the images of the electrons (each line is characterized by the Coulomb 
energy l/4d) and by solid lines the attractive interactions between the electrons and their 
images (each with the energy —l/4d). Now, in the case of the D~ , due to the presence of 
the two electrons (see Fig. [1]), we have twice the energy shift of the single electron problem 
2AE = {l/4d), and there is also twice a negative shift due to the interaction of each electron 
with the image of the other electron 2{—l/A)d (see Eq. (I5bl) ). As a result these energy shifts 
compensate each other, and we find that the energy for d > 10a b is equal to the energy of 
the D~ without image charges and is very close to the energy Ehuik/2R* = —0.5259 of D~ 



in the bulk [l|]. The same compensation takes place in the case of a semiconductor-dielectric 
interface where the large d contribution of image terms in the energy of the neutral D° center 
is about —Q/M 



The average position of the electrons in the z-direction and their extend in the plane 
parallel to the interface are shown in Figs. [5]^a) and[5]^b), respectively. Notice that one of 
the electrons, also called the inner electron, follows very closely the behavior of the electron 
bound in the neutral D^ system. This is very similar to what was found previously for the 
case of a bulk D~ 7|]. The second electron, called the outer electron, is more extended in 
the p-plane, i.e. its average value (p2) is about three times larger than for the inner electron. 
Notice that (p2) for the outer electron increases rapidly with decreasing d when d < 2aB 
signaling a rapid decrease of the binding energy and ultimately an unbinding of the outer 
electron. For d > 4.750^ the outer electron is attracted towards the interface, i.e. (2:2) < d, 
as a consequence of the image charges which is responsible for the enhanced D~ binding 
energy. For smaller d values we have (^2) > d and the outer electron is pushed away from the 
interface mostly as a consequence of the boundary condition at the interface. This behavior 
is also illustrated in Fig. [6] where we show the contour plots of the electron density of the 
outer and inner electron, i.e. |'?/'2(^2,P2)P = p(-2i)V(-22)^5'(?i, ?2)^V^(ri, ^2)^ with ri, pi and 
{61—62) taken as their average value, for d/as = 2, 5, 9. Notice that the electron distribution 
is asymmetric when the donor is close to the interface and that in such a case a large part of 
the distribution is found with z > d. In order to illustrate the effect of the electron-electron 
correlation we show in Fig. [7] the conditional probability P(r, ro) = |'?/'(r', ro)| • This is the 
probability to find an electron at position r when the other electron is fixed at position 
tq. We fix one of the electrons close to the donor (i.e. it is the inner electron) and put 
it in three different positions with respect to the interface-donor axis. Notice that the 
electron: 1) has the highest probability to be close to the donor, 2) it is repelled by the fixed 
electron, and 3) it has a non-zero probability to be located at tq. The reason is that the 
two electrons have opposite spin and therefore the Pauli exclusion principle is not applicable. 

B. Semiconductor-dielectric interface 

In this subsection we investigate the ground state energy of a D~ center near a 
semiconductor-dielectric interface. We used material constants for the semiconductor side 
corresponding to 63 = 11.9 (Si) and for the oxide side eq = 3.4 (Si02). 

The potential energy between the particles is still given by Eq. ([1]) with corresponding 
value oi Q = 0.556. The results are presented in Fig. [8] for the ground state energy and in 
Fig. inifor the binding energy. The second electron is bound to the D^ for d > S.das which 
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compares to d = 3.2aB when we ignore the image charges. This can be explained by the fact 
that for the D~ system near a semiconductor-dielectric interface the two electrons induces 
two negative image charges while they are only attracted by one positive image charge 
coming from the donor. In Fig. [S] we show also the energy of the D^ by the dashed blue 
curve which corresponds to the situation where the second electron is at infinity. Notice 
that both curves cross at d = 6.80^ and for smaller (i-values the D~ state has a higher 
energy and consequently the second electron will be unbound. The binding energy of the 
second electron to the D° is plotted in Fig [9] by the solid red curve. This result approaches 
the bulk result slowly and from below. The inset of Fig. M shows the large d-behavior which 
is fitted by the curve Eh/2R* = 0.0259 — OAAdas/d. In Fig. dUlwe show the electron density 
of the inner and outer electron at d/as = 5, 6.3 and 8, corresponding to the black squares 
marked in Fig. |8l Notice that the outer electron is much further extended in space than the 
inner electron and is pushed away from the interface. 

C. Semiconductor- vacuum interface 

For the D~ center near a semiconductor-vacuum interface, the material constant is chosen 
to be Es = 11.9 for the semiconductor side and Eq = ^ for the vacuum side, leading to 
Q = 0.845. 

The qualitative behavior for the energy of the D~ and D° as function of d is similar 
to the one shown if Fig. [HI Due to the larger Q, the role of the image charges becomes 
more important and the outer electron becomes bound for d = 4a b- The crossing point 
between the D~ and D^ curves is pushed to d = 8.80^. The binding energy is shown in 
Fig. [9] by the dashed blue curve. Notice the strongly reduced binding energy as compared to 
the Si/Si02 interface and even more so when we compare it with the semiconductor-metal 
interface case. The large d-behavior is shown in the inset which could be fitted by the curve 
Eb/2R* = 0.0259 - 0.2l8aB/d. 

IV. CONCLUSION 

We proposed a variational approach to investigate the energetics and the wave function 
extend of the spin singlet ground state of the D'^ system that is located near a semiconductor- 
metal or a semiconductor-dielectric (vacuum) interface. As a trial function we used a mod- 
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ified Chandrasekhar type wave function, which differs with the Chandrasekar variational 
two-electron wavefunction in the following way: 1) it satisfy the boundary condition at the 
interface and 2) it takes into account all the Coulomb interactions with the image charges. 
This makes the wavefunction no longer spherical symmetric. This variational approach gives 
the well-known rf — )> oo limit. 

We obtained a nonmonotonic behavior of the D~ binding energy as function of donor 
position near a semiconductor- metal interface with a local maximum for a donor distance 
from the interface of about S.Sa^. For smaller rf-values the binding energy decreases and 
for d < l.las the D~ becomes unbound which is mostly a consequence of the bound- 
ary condition of the two electron wavefunction at the interface. At large d the outer 
electron is attracted to the interface because of the positive total image charge. For a 
semiconductor-dielectric (vacuum) interface the D~ binding energy is strongly reduced and 
is a uniform increasing function of d. When the donor moves towards the interface, the 
energy of the D~ reduces and the system becomes unbound for d < S.Sa^ {d < Aas) near a 
Si/Si02 (Si/vacuum) interface. While for a neutral D^ center near a semiconductor-metal 
(semiconductor-dielectric) interface the contribution of the image terms to the binding en- 
ergy is approximately given by the expression AE = l/4d {—Q/Ad) for large d. In the case 
of the D^ system, due to the presence of the two electrons, a complete compensation of such 
terms takes place. 

The present calculation was done within the effective mass approximation which is ex- 
pected to be valid for the considered length scale d > as- For a donor very close to the 
interface, i.e. within three monolayers, the deformation of the lattice close to the interface, 
i.e. strain relaxation near the donor, may invalidate the effective mass approximation. We 
also neglected the penetration of the electrons in the metal (and the dielectric) which for the 
obtained binding energies is expected to be a good approximation. In order to go beyond 
the present effective mass approximation one can use approaches such as the tight-binding 
approximation or approaches based on the density functional theory. 
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FIG. 1: (Color online) Schematic picture of the interaction between the D~ system and its image 
charges for the case of a semiconductor-metal interface. Dashed (solid) lines indicate the repulsive 
(attractive) interactions between the electrons, the donor and their images. 
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(a) semiconductor-metal interface 

for the first electron 

for the second electron 




(b) semiconductor-dielectric interface 
for the first electron 
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FIG. 2: (Color online) The mean-field potential energy along the z-axis (with p = 0) for a donor 
located at d/as = 5 in case of (a) a semiconductor-metal and (b) a semiconductor-dielectric 
{es = 11.9, Ed = 3.4) interface. We show the result for the potential of a single electron (dashed 
blue curve) and for the electrostatic potential due to the unperturbed D^ system with the electron 
in its ground state (solid red curve). 
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FIG. 3: (Color online) The ground state energy (solid red curve) for a D~ system near a 
semiconductor-metal interface vs. the distance of the donor d/uB from the interface. Dashed 
green curve is the results when the image charges are neglected and by black stars we indicate 
the results when an additional variational parameter 6z is included in the e-e correlation term (see 
Eq. dS])). For comparison the energies for a neutral donor (D^) (dot-dashed blue curve) and in the 
absence of the image charges (dot-dot-dashed dark red curve) are also shown. 
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FIG. 4: (Color online) The binding energy (solid red curve) for a D~ center near a semiconductor- 
metal interface vs. the donor- interface distance djas- The binding energy in the absence of S/ (the 
energy of an electron bound to the interface) is shown by the dashed blue curve. The horizontal 
dashed black line shows the bulk limit of the binding energy. The inset shows the behavior of the 
binding energy for large d (black circles) which is fitted to the curve E[)/2R* = 0.0259 + 0.255aB/d 
(dashed blue curve). 
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FIG. 5: (Color online) For the semiconductor-metal interface: (a) The average distance between 
the electrons and the donor along the z direction vs. the donor position d. (b) The average value 
of the electron positions along the p direction vs. the donor position d. The inset of figure (a) 
shows the behavior of (zj — d) at large d. 
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FIG. 6: (Color online) The electron probability density of the D~ center near a semiconductor- 
metal interface for different donor distances d/uB = 2, 5 and 9 when the other electron is fixed to 
its average position (see the corresponding average values in Fig. [5]) as well as the angle between 
two electrons (^i — 62) is set to its average value. The upper three figures correspond to the inner 
electron density (when the outer electron is fixed) , whereas the lower three figures are densities of 
the outer electron. Blue (red) area represents low (high) probabilities and the white area represents 
almost zero probability. 
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FIG. 7: (Color online) The conditional electron probability density of the outer electron where the 
inner electron is fixed when the D~ center near a semiconductor-metal interface with d/as = 3. 
The white circles marked out the location of the inner electron, that is (from left to right, in the 
form (x/as, y/as, z/aB)) (0, 0, 2), (1, 0, 3) and (0, 0, 4). 
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FIG. 8: (Color online) The ground state energy (solid red curve) of a D~ center near a 
semiconductor-dielectric interface vs. the donor-interface distance d. The energy for a neutral 
donor (-D") (dashed blue curve) near a semiconductor-dielectric interface is shown for comparison. 
The black squares are the d- values for which we show the electron density in Fig. [101 



22 




d/ap 



FIG. 9: (Color online) The binding energy for a D~ center near a semiconductor-dielectric 
(semiconductor-vacuum) interface vs. the donor distance d/as indicated by solid red (dashed 
blue) curve. The horizontal dashed black line shows the bulk limit of the binding energy. The 
inset shows the behavior of the binding energy (symbols) for large d. The solid red (dashed blue) 
curve is a fitting curve based on the equation E/2R* = 0.0259 — xas/d with the fitting parameter 
X = 0.145 (x = 0.218). 
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FIG. 10: (Color online) The electron density of the D~ center near a semiconductor-dielectric 
interface at donor distances d/as = 5, 6.3 and 8 (see the black squares in Fig. [8]). The top three 
figures correspond to the inner electron density (when the outer electron is fixed to its average 
position), whereas the bottom three figures are the densities of the outer electron with the inner 
electron fixed to its average. Blue (red) area represents low (high) probabilities and the white area 
represents almost zero probability. 
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FIG. 11: (Color online) Similar as Fig. [8] but now for the D and D^ near a semiconductor- vacuum 
interface. 
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